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ABSTRACT 

The desorption characteristics of molecules on interstellar dust grains are important 
for modelling the behaviour of molecules in icy mantles and, critically, in describing 
the solid-gas interface. In this study, a series of laboratory experiments exploring the 
desorption of three small molecules from three astrophysically relevant surfaces are 
presented. The desorption of CO, O2 and CO2 at both sub-monolayer and multilayer 
coverages was investigated from non-porous water, crystalline water and silicate sur- 
faces. Experimental data was modelled using the Polanyi-Wigner equation to produce 
a mathematical description of the desorption of each molecular species from each type 
of surface, uniquely describing both the monolayer and multilayer desorption in a sin- 
gle combined model. The implications of desorption behaviour over astrophysically 
relevant timescales are discussed. 
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INTRODUCTION 



The interaction of molecular species with surfaces is of crit- 
, ical importance in astrophysical environments. Much in- 
terstellar chemistry occurs on or in the icy layers which 
cover dust grains in molecular cloud s, the birthplace of 
, stars and planets dFraser et al.l [2002). These dust grains 
are believ ed to be com posed of silicates and carbonaceous 
'material (|Drainei 12003). The silicate is at least ~ 95 % 
amor phous, as determin ed from the breadth of the 9.7 /im 
band dLi fc Drainei r2002), and although its composition has 
not been determined precisel y, is well approximated by 
amor phous olivine (MgFeSi04) IjSofia fc Mey er 2001; DraineJ 
2003). The silicates are believed to be 'fluffy' in natu re, with 
a larg e surface area upon which reactions can occur l|Mathis! 
1 19981 ). 

At the low temperatures in molecular clouds (T < 20 K) 
atoms and molecules freeze out onto these silicate surfaces, 
forming the icy mantles. Amorphous solid water (ASW) is 
by far the largest component of this ice, with abundances 
of ~ 1 x 10~ 4 with resp ect to the total H column density 
(Williams & Hcrbst 2002), equivalent to coverages of up to 
100 monolayers (ML). Given that t he extinction thresh old 
for H 2 mantles is Ay ~ 3.3 mag |Whittet et al.lll98Sl ), it 



E-mail: francois.dulieu@obspm.fr 



is reasonable to assume that between the dense pre-stellar 
cores (where dust grains are completely coated by this icy 
mantle) and the cloud edges (where competition between ice 
formation and photod esorption of H2O y ields a population 
of bare silicate grains (jSmith et al.lll993l )). there must be a 
region where icy surfaces and bare silicates co-exist. 

As interstellar regions evolve, the icy grains are fur- 
ther processed, eit her by gentle heati ng or cyclic desorption- 
deposition events IjVisser et al.ll2009l ). producing crystalline 
H2O. This has been detected in various stellar and pre- 
planet ary environments, i ncluding, for example, M giant 
stars (|Omont et al.l Il990l ). Quaoar in the Kuiper Belt 



Jjewitt fc Luull2004l). Trans Nept une Objects (|Merlin et all 
I2007D . comets dLisse et al.l 1200a), and in th e outer disks 
around T Tauri stars (|Schegerer fc Wolj|2010l ). 



As gas-grain modelling of interstellar environ- 
ments becomes more sophisticated ([Wakelam et alj l2010l : 
lAcharvva et~aH 1201 ll : ICuppen fc Garrodl l201ll )~ key ques- 
tions impeding the full implementation of surface chemistry 
include: what effect does the underlying grain surface have 
on the desorption characteristics of key molecules, and to 
what extent is this desorption affected as we move from the 
multilayer to sub-monolayer coverage regime? A third issue 
is how to realise the transition between the multilayer and 
monolayer regimes in a gas-grain model without overloading 
the model with many more layers of complexity. 



© 2002 RAS 



2 J. A. Noble et al. 



CO, O2 and CO2 are all interstellar molecules that 
potentially could populate bare interstellar grains at sub- 
monolayer coverages. On ice, only one previous study fo- 
cussed on sub-monolayer coverages of CO, highlighting the 
spectroscopic rather than desorption charact eristics of the 
porous ASW:CO system l|Collings et alj|2005h . To date, no 
desorption studies of CO, O2 or CO2 have been made on a 
silicate surface, at either multilayer or sub-monolayer cover- 
ages. 

CO is th e second most abundant interstellar molecule 
l|Tielenj|2005h , and is known to form in the gas phase, then 
freeze-out on to H2O covered grains to form overlayers of 
pure CO ice, the layer thickness criticall y depending on gas 
densi ty rather than grain temperature (|Pontoppidan et al.l 
200j|) • Consequently, extensive temperature programmed 
desorption (TPD) studies have been made of multilayer CO 
coverages on various surfaces, including Au llAcharvva et al 



20071 ; iBisschop et al.ll200rj; ICollings et & . I l2004l ; I Fuchs et al 



20061). H2O (ICollings et al.l l2003al lbl), a meteorite sample 
( Mautner et al-luOOq) and high ly oriented pyrolytic graphite 
(HOPG) (jUlbricht et al.|[2006l ). 

It is widely accepted that both CO2 and H2O form 
on dust grains in molecular clouds, and recent stud- 
ies suggest that at le ast some of the CO2 and H2O is 
form e d concurrently ([Goumans et al. I 120081 ; I Noble et al.l 
l201ll ; lloppolo et al.ll201ll ). Since CO is the key precursor to 
CO2 formation, such mechanisms would require CO freeze- 
out onto bare grain surfaces long before CO ice, or even large 
quantities of H2O ice, are detected. Although tenuous, spec- 
troscopic observational a nd experimental ev idence exists for 
such a freeze-out process l|Fraser et al.ll2005l ). but an investi- 
gation of the desorption behaviour of CO at sub-monolayer 
coverages on ice and silicate surfaces is vital if we want to 
be certain that it can reside at an interstellar grain surface 
long enough to form more complex species. 

Similarly, O2 is a potential precu rsor in H2O formation 
jloppolo et al.ll200Sl ; lOba et all 120091 ), and is also a species 
likely to form in the gas phase then freeze-out onto grain 
surfaces, rather than forming on the grain itself. However, 
as a homonuclear diatomic O2 is infrared inactive, so its de- 
tection in interstellar i ce, though occasio nally claimed via 
a forbidden transition (|Elsila et al.l Il997h remains elusive. 
It is not clear whether this is related to the weak transi- 
tion probability, lack of a significant O2 population in the 
ice, or that the O2 has rapidly reacted (upon adsorption) 
to exclusively form H2O. Nevertheless, the multilayer des- 
orption behavio ur of O2 has been studied previously alo ng- 
side CO and N 2 l|Acharvva et al.ll2007l;lFuchs et al.ll2006h . as 
well as on Au, porous ASW jCollings et al.ll20ol ) and Ti0 2 
jDohnalek et al.ll2006h . These studies show that the multi- 
layer desorption characteristics of O2 are very similar to CO, 
so it is interesting from both a chemical and astrophysical 
viewpoint to also investigate the desorption behaviour of O2 
at sub-monolayer coverages on ice and silicate surfaces. 

The desorption characteristics of CO2 have not been 
extensively studied on any surface, despite it being one of 
the most abundant solid phase molecular species in the in- 
terstellar medium. Recently the desorption characteristics of 
multilayer CO2 were reported from p orous ASW, and hav e 
been previously studi ed on HOPG (Tuibricht et al.l l2006h . 

However, given 
and experimen- 



porous ASW and Au (Collings et al 
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tal jNoble et alj|201 ll ) evidence shows that a fraction of CO2 
ice must form concurrently with the water ice layer, some 
CO2 molecules must populate both the bare silicate grains 
and the ice layers at sub-monolayer coverages. 

Here we present an experimental study of the desorption 
of CO, O2 and CO2 from three different surfaces; non-porous 
ASW, crystalline ice, and amorphous olivine-type silicate. 
As these are, for the first time, all u ndertaken in the sam e 
experimental set-up - FORMOLISM l|Amiaud et al.ll2006l ) - 
we are able to investigate both the individual effect of each 
surface on the desorption characteristics, as well as deter- 
mining whether the molecular composition or morphology 
of the surface is most relevant in determining the desorp- 
tion behaviour of molecules. For each of the nine combina- 
tions of CO, O2 and CO2 on each surface, our study has en- 
compassed both the multilayer and sub-monolayer regimes. 
These data are modelled to determine a simple analytical ex- 
pression which accurately calculates both the sub-monolayer 
and the multilayer desorption energies. By changing the 
model to incorporate interstellar, rather than experimental, 
heating rates, we are able to address the key questions above, 
namely: what effect does the underlying surface have on the 
desorption characteristics of the molecules adsorbed there, 
and how is this desorption modified in the sub-monolayer 
coverage regime? 



2 EXPERIMENTAL 

The experiments were conducted using the FORMOLISM 
set-up (FORmation of MOLecules in the InterStellar 
Medium), described in detail elsewhere (|Amiaud et al 1 120061 : 
lLattelais et al.ll201ll ) . Briefly, the set-up consists of an ultra- 
high vacuum chamber (base pressure ~ 10~ 10 mbar), con- 
taining a silicate-coated copper sample surface, operating at 
temperatures between 18 and 400 K. The system is equipped 
with a quadrupole mass spectrometer (QMS), which is used 
for the TPD experiments. A sample of either O2, CO or 
CO2 was deposited onto the surface at 18 K via the triply 
differentially pumped beam line; a linear temperature ramp 
was then applied to the surface, and the QMS used to mea- 
sure the desorption of each species into the gas phase, as a 
function of temperature. 

Three surfaces were investigated in this study: non- 
porous ASW (H20( np j), crystalline ice (H20( c j) and amor- 
phous silicate (SiO^). The SiO^ surface w as recently in- 
stalle d in the FORMOLISM experiment |Lemaire et al.l 
l2010h . and mimics bare dust grains in molecular clouds. 
The silicate is amorphous in nature, as evidenced by in- 
frared spectroscopic studies, while TPD experiments, in- 
cluding those presented in this study, reveal the surface to 
be non-porous on the molecular scale. 

The silicate surface is capable of reaching base tem- 
peratures of 18 K, which is how the molecular deposition 
temperature on all the surfaces was pre-determined. For the 
ice surfaces, 50 ML films were grown on top of the silicate 
surface by spraying water vapour from a microchannel array 
doser located 2 cm in front of the surface. The water vapour 
was obtained from deionised water which had been purified 
by several freeze-pump-thaw cycles, carried out under vac- 
uum. H20( np ) mimics the ASW which comprises the bulk of 
interstellar ice, and H20( c ) mimics the crystalline ice seen in 



© 2002 RAS, MNRAS 000,riHT3l 



Thermal desorption characteristics of CO, 2 and C0 2 on H 2 0( np ), H 2 0^ and SiO x . 3 



Table 1. Description of experimental exposures. 



Molecule 


Surface 


Depositions" 
ML 


o 2 


H 2 (c) 
SiO* 


0.25, 0.50, 0.75, 1, 1.25 c , 1.38, 1.75, 
2.25 c , 2.75 

0.08, 0.17, 0.50, 0.67, 1, 1.33 
0.20, 0.60, 1, 1.20, 1.60, 1.67 


CO 


H20( np ) 

H 2 {c) 
SiO* 


0.25, 0.50, 0.75, 0.88, 1, 1.06\ 1.13 6 , 
1.25 6 , 1.75 6 , 2.50 

0.10, 0.20, 0.60, 0.80, 1, 1.20 6 , 1.60 
0.50, 1, 1.25 6 , 1.75, 2.25 


C0 2 


H20( np ) 

H 2 (c) 

SiO* 


1, 2 b , i b , 6 b , 10 b , 20 
1, 2, 3, 5 C , 10, 20 
0.5, 1, 2, 3, 5, 10 



a defined as e ^^ e , where N mono is defined in $ |3~T1 
6 Omitted from model (shown in grey in Figure [T}. 
c Modelled, but omitted from Figure \3\ for clarity. 



some star-forming regions. To produce H20( np ), water was 
dosed while the surface was held at a constant temperature 
of 120 K. To form H20( c ), the surface was held at 120 K 
during the deposition, then flash heated at 50 Kmin" 1 to 
140 K, and finally at 10 Kmin" 1 to 142.5 K. For each type 
of ice surface, the temperature was then held constant until 
the background pressure in the chamber stabilised, before 
cooling it back down to 18 K, at which temperature adsor- 
bates were dosed onto the respective surfaces. 

For each of the nine combinations of molecule and sur- 
face, a series of depositions were made, varying the adsor- 
bate surface coverage, from sub-monolayer to multilayer, as 
outlined in Table [T] After the first adsorbate exposure, the 
surface was then heated at a rate of 10 Kmin -1 until the 
adsorbate was fully desorbed from the surface (below 100 K 
on the ice surfaces and 130 K on the silicate surface). The 
original surface therefore remained intact, and by purpose- 
fully keeping the temperature of the ice below 100 K no 
further thermal annealing of the ice occurred during an ex- 
perimental sequence. Subsequently, the surface was then re- 
cooled to 18 K, and a new coverage of the same adsorbate 
added. In this way we could be certain the same surface 
was used for a single series of adsorbate exposures, allow- 
ing real comparison of the coverage effect for the first time. 
In the multilayer, the underlying surface is generally not 
important to the desorption, as adsorbate-adsorbate inter- 
actions between molecules determine the kinetics. In the 
sub-monolayer, the interaction between the adsorbate and 
the substrate is fundamental in determining the desorption 
characteristics. Thus, by including both regimes, we hoped 
to probe the differences between these two coverages for dif- 
ferent molecules and different surfaces. 

A porous ASW surface was not included in this study 
for two reasons. First, the presence of pores in the ice sur- 
face can produce very complex desorption profiles, due to 
the competing kinetics w hich arise after entra pment of ad- 
sorbates in the pores (e.g. lCollings et alj|2003ah . Worse still, 
pore collapse starts at only 10-15 K, so the experimental 
approach used here, of exposing a single ice surface to a va- 
riety of adsorbate coverages, is simply not possible with the 
porous ASW. A new ASW sample would have been required 



for each experiment; each time we form a porous ice in the 
laboratory its structure will be slightly different, resulting in 
a massive change in surface area and non-reproducible ad- 
sorbate exposures being required to obtain monolayer cover- 
ages. Furthermore, one aim of this experiment is to compare 
surface characteristics using a single model for the whole 
dataset, and such an approach would not have been possi- 
ble when including both porous and non-porous ices, again 
because the monolayer surface area is difficult to define. Sec- 
ondly, it has recently been shown experimentally that it is 
very likely that al l water ice in inters tellar regions will be 
non-porous ASW (|Accolla et al.ll2~01ll ). even though obser- 
vations previously suggested that invoking a porous ASW 
was the only way to explain observed qua ntities of mixed 
H2O and CO ices (jPontoppidan et al.l l2003). 



3 RESULTS AND DISCUSSION 
3.1 Experimental data 

All the TPD results are shown in Figure [1] The results are 
ordered top to bottom by surface type, H20( np ), H20( c ) 
and SiOz, and left to right by adsorbate molecule, O2, CO 
and CO2. For each molecule-surface combination, the main 
window shows the series of sub-monolayer exposure TPD 
spectra, derived from the lowest molecular coverage, up to 
1 ML (as per the values given in Table [IJ. The inset shows 
the equivalent data for coverages of 1 ML and above. 

The spectrum which represents 1 ML is chosen by visual 
inspection of all TPD spectra in the series, and is defined 
in this study as the highest coverage spectrum which in- 
cludes only (sub-)monolayer desorption characteristics, and 
has no multilayer component. For a fully wetting molecule, 
such as CO or O2, this coverage is assumed to contain 
10 15 molecul es cm" 2 , the genera lly accepted definition of 
a monolayer l|Amiaud et al . 2006). In the case of CO2 , how- 
ever, a lower exposure of the species was required before 
intermediate and multilayer desorption characteristics were 
observed i.e. the 1 ML coverage contains fewer than 10 15 
molecules. Thus, in this study, we have defined the cov- 
erage as 6 = e ^P oaure ; where exposure is the number of 
molecules deposited on the surface and N mono is the number 
of molecules in the designated 1 ML deposition. This defini- 
tion is necessary to allow modelling of the CO2 data in the 
same fashion as the CO and O2 data, as discussed in § 13.21 

By inspection of Figure [1] it is immediately apparent 
that, for all cases, the TPD line profile differs between the 
sub-monolayer and the multilayer regime, as expected. In 
some cases (shown in grey on Figure [TJ, an intermediate re- 
gion is visible in those TPD spectra which trace coverages 
between 1 ML and the multilayer, for example in CO or CO2 
on H20( n p). It is clear that structural reorganisations must 
occur at these intermediate coverages, as desorption spectra 
with double and triple peaked profiles are observed. Describ- 
ing such processes is beyond the scope of this paper and for 
this reason these particular data have been omitted from the 
later modelling, only being shown here for completeness. 

For both O2 and CO2 it is also apparent from Figure [1] 
that the TPD line profiles are similar in the sub-monolayer 
regime, independent of the underlying surface. For CO this 
is clearly not the case, although the multilayer data are 
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Figure 1. Temperature-programmed desorption spectra of O2, CO and CO2 from H20( nJ) ), H20( c ) and SiOa;. For each molecule-surface 
combination, the main window shows a series of TPD spectra from the lowest sub-monolayer coverage up to 1 ML; inset are TPD spectra 
from 1 ML up to the highest multilayer coverage, as specified in Table[T] Spectra shown in black were included in the model, whilst those 
in grey were omitted. In all cases, the line profile of the TPD spectra in the multilayer regime differ from those in the sub-monolayer 
regime and, with the exception of CO, the form of the TPD profiles is very similar in the sub-monolayer regime, independent of surface. 
The desorption profiles of all three molecules from the Si02 surface in the sub-monolayer regime strongly resemble each other, unlike 
the ices. See text for full discussion. 



more consistent across all surfaces. Likewise, in both the 
sub-monolayer and multilayer regimes, the TPD spectra for 
different adsorbates on the SiO^ surface strongly resemble 
each other; the same can not be said for the ice surfaces. For 
all the adsorbates, it is evident that multilayer desorption 
from the crystalline ice surface starts at a slightly higher 
temperatures than from either the amorphous ice or amor- 
phous Si02 surface, suggesting a trend dependent on the 
morphology of the underlying surface. Specifically, for O2, 
the peak temperature occurs between 30.9 and 32.6 K, but 
starts around 24-25 K on amorphous surfaces (H20( np ) and 
SiOrr), compared to 25-26 K from H20( c ). For CO, the peak 
temperatures are in the range 28.1-28.9 K, desorption starts 
at 22-23 K from the amorphous substrates and 23-24 K 
from H20( c ). For CO2, the peak desorption range is wider 



(79.3-83.0 K) and desorption begins at 63-65 K from SiO^, 
64-65 K from H 2 (np ) and 66-67 K from H 2 (c) . 

To compare the sticking efficiency of each molecule on 
each surface, th e area under the T PD peak, corrected for 
QMS sensitivity l|Matar et al.ll20l"(il ). was plotted against de- 
position time for each molecule on each surface, see Figure[2] 
Data can be divided into adsorbates - circle (CO), triangle 
(O2), and square (CO2) - and surfaces - H20( np ) in black, 
H.20( c ) in light grey and SiO^ in dark grey. Such a direct 
comparison is possible because the experimental method em- 
ploys the same surface for a series of adsorbate exposures. 
Furthermore, the dosing rates of all three molecules, and 
the effective beam pressures employed, are always approx- 
imately equal. Straight lines of the form y = mi + c were 
fitted to each molecule-surface combination and overplot- 
ted on the experimental data in Figure [2] Within experi- 
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Figure 2. Area under TPD peak, corrected for QMS sensitivity, 
plotted against deposition time. Experimental data is plotted as 
points, with H20( np ) in black, H20( c j in light grey and SiO^ in 
dark grey; CO as circles, CO2 as squares and O2 as triangles. 
Overplotted are the straight line fits to the data, as described 
in the text. Colours are as above for the three surfaces, with 
CO solid, CO2 dashed and O2 dot-dashed lines. As expected, 
for each molecule-surface combination, the area increases with 
deposition time, and the straight line fit goes through zero, within 
experimental errors. For each molecular species, sticking occurs 
approximately cquivalently on all three surfaces, with an error 
of < 10 %, although sticking is slightly lower on H20( np ) for all 
molecules. 



mental uncertainty (shown as error bars in Figure |2j, these 
straight line fits all pass through zero, as expected, indicat- 
ing that the sticking probability is constant (within < 10%) 
across the dynamic range of each experiment. This means 
that no discernible change in sticking probability occurs for 
any molecule on any surface between the sub-monolayer and 
multilayer regimes. Adsorbate coverage is always linear as a 
function of exposure time. Crucially, for all three adsorbate 
molecules, the sticking efficiency on H2 0( np ) (black lines in 
Figure [2]) is slightly lower when compared with the other 
two surfaces. This is likely due to fact that the average ad- 
sorption energy on H20f np ) peaks at a lower energy than for 
the other substrates. 

One final point to emphasise is the distinction between 
the time at which a surface will have been exposed to 
an equivalent dose of 10 15 molecules (generally defined as 
the number of molecules in one monolayer for a wetting 
molecule), and the coverage (exposure time) at which the 
TPD data show a deviation from sub-monolayer to mul- 
tilayer behaviour, defined above as N morlo . For CO and 
O2, Nmono is reached on all surfaces between around 250- 
300 s exposure; for CO2 the switch from monolayer to multi- 
layer occurs at much lower surface coverages (exposure times 
~ 30 % those for CO and O2). We postulate that this dif- 
ference is related to the wetting behaviour of the molecules, 
and come back to this issue in S 13.31 



3.2 Modelling 

To obtain a more quantitative understanding of the interac- 
tions between the molecules and the surfaces, the TPD spec- 
tra were reproduced from an empirical kinetic model. To do 
this, the results were divided into sets of TPD spectra as- 
sociated with sub-monolayer and multilayer behaviour and 
the kinetic parameters in each regime determined indepen- 
dently, before being recombined to reproduce the experimen- 
tal TPD data. A simple model was therefore derived to fit 
all nine molecule-surface combinations from sub-monolayer 
to multilayer adsorbate coverages. 

On any surface, at any coverage, thermal desorption can 
be described in terms of an Arrhenius law, expressed by the 
Polanyi-Wigner equation, where the rate of desorption, r, is 
described by: 



dN 



= AN"e 



-E ad ./kT 



(X) 



where A is the pre-exponential factor, N is the number 
of adsorbed molecules on the surface (cm -2 ), n is the or- 
der of the reaction, E a ds is the energy of adsorption of a 
molecule to the surface (eV), k is the Boltzmann constant 
(eV molec _1 K~ 1 ), and T is the temperature of the surface 
(K). The units of A depend on n: molecules 1 ~ n cm~ 2+2n s~ 1 . 
Reformulating Equation [T] to reflect the measured experi- 
mental TPD signal gives: 



dN_ A 

dT p 



,/kT 



(2) 



where p is the rate of heating, in this case dT/dt 
10 Kmin" 1 . 

When the rate of desorption is independent of surface 
coverage, as is the case in multilayer desorption, n — and 
the TPD desorption profiles share a common leading edge 
as is readily apparent from the insets in Figure [T] 

For each molecule-surface combination, the leading edge 
of the highest multilayer exposure TPD spectrum was used 
to determine E a ds (in K). Although some multilayer des- 
orption TPD data show tha t the order of reaction can de- 
viate slightly from zero (e.g. iBrown fc Bolinal l|2007l )) it is a 
reasonable approximation, so henceforth w e assu me n — 0. 
As applied previously by lAcharvva et al.l (|200?) , the pre- 
exponential factor in Equation [5] was assumed to be a func- 
tion of E a ds approximated by: 



Nml ■ v = N M l 



2NMLEgd S 

7T 2 M 



(3) 



where M is the mass of the adsorbate molecule, and 
Nml ~ 10 15 cm -2 . The key advantage of this method 
is that the multilayer fitting requires only one variable, 
Eads, removing the interdependency of fitting both the pre- 
exponential factor and E a ds concurrently. The results of fit- 
ting -Eads are presented in Table [2] 

In the sub-monolayer regime, the rate of desorption is 
dependent on surface coverage, 6 = N/N morlo ^ 1, where 
N morlo (as defined in § I3.1|l is the maximum number of 
molecules on the surface prior to the onset of intermedi- 
ate or multilayer desorption behaviour. Previous publica- 
tions have s hown that the reaction o rder can deviate slightly 
from 1 (e.g. IBrown fc Bolinal (|2007f )1. but for the purposes 
of this work we assume n = 1. In general, in the case of 
first order desorption, the peaks of the TPD data appear 
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Table 2. Multilayer desorption parameters calculated using 
Equation [2] compared to previous literature values in italics 



Molecule 


Surface 


Calculated values 








V 


Eads 






10 26 molec cm _2 s _1 


K 




H20( np ) 


6.9 


898(30) 


o 2 




i X) 






SiOa; 


6.9 


895(36) 




Au 


6". 9 


912(15) 








Oil (r /(H ;r 1c 

y^ci ( ^0 y 






7.1 


828(28) 




H 2 (c) 


7.1 


849(55) 


CO 


SiO* 


7.1 


831(40) 


Au 


7.2 


858(15) b 
826(24) d 








855(25)" 






9.3 


2267(71) 




H 2 (c) 


9.5 


2356(83) 


co 2 


SiO* 


9.3 


2269(80) 




porous ASW 




2690(50) f 




HOPG 




2982(-)b 



a For values calculated in this work, the values in parentheses 
are the total calculated error, including the 3 a statistical 
errors on E aaa calculated during the fit of Equation [2] and the 
experimental uncertainties. In all other cases the error is as 
q uoted in the relevant l iterature. 

b lAcharwa et al ] j2007l ). using an identical multilayer modelling 
method. 

c iFuchs et alj l|2006h , calculated assuming a fixed pre-exponential 
factor. 

^Collings et al ] j2003bh . calculated assuming a fixed pre- 
exponential fa ctor. 

Hdberg et all j2005h : IFuchs et all [|2006h : iBisschop et~all I^OOrl ). 
c alculated assuming a fixed pre-e xponential factor. 
■^ ISandford &: All amandola ( 1990), calculated using spectroscopic 
d ata, assuming a z e roth o rder desorption profile, 
f iBurke fc Brownl (|201(J) . calculated empirically with a non- 
integer reaction order. 



at a single temperature value, pro vided that the des orption 
energy is independent of coverage (Fras e~r~et al.ll20o"lh . How- 
ever, as is evident from the main panels in Figure [1] in this 
instance it is the trailing edge of the TPD curves that are 
common in each molecule-surface system. Such behaviour 
is known to arise when the adsorbate is occupying multi- 
ple sites on the surface, such that the adsorption energy is 
a continuou s function of the number of molecules on the 
surface (e.g. iDohnalek et al.ll200ll ; lAmiaud et all 120061 ). By 
inverting Equation [2] a TPD spectrum can be converted to 
a function, E(N), given by: 

E(N) = -kTln^L. (4) 

To fit E(N) for each molecule-surface system, we first 
select the TPD spectrum equating to the surface exposure 
just prior to the appearance of multilayer peaks, or (in the 
case of CO2) more complex desorption characteristics, and 
define this as the monolayer coverage, i.e. — N/N mono = 1, 
ensuring the maximum range of adsorption sites and ener- 
gies at the surface are included. This is a very reasonable 
approach, given that it is well known that E(N) is indepen- 



dent of reaction order at exactly N = N mono (lUlbricht et alj 
2006). This monolayer TPD is then fitted using an analytical 
expression of E(N), redefined in terms of surface coverage, 
E(6), given by: 

E{9) = Q + Qi 9 + OL 2 e 2 + a 3 e" 4 e + 

(5) 

Q5 e + c*7 e B + ag e . 

The resulting coefficients ao - aio are listed in Table [3] 
and can then be used to calculate E{6) from Equation [5] 
which can then be substituted into Equation [3] along with 
an appropriate value for A, to reproduce all TPD data for 
a specific molecule-surface system. 

As the choice of A has a large effect on the value of E(8) 
calculated from the model, we first attempted to reproduce 
the sub-monolayer cov erages by concurrently optimising A 
Iptirniman et al.lll996T ). During the calculation of E(N), A 
was set at a fixed range of values from 10 5 - 10 15 s _1 . The 
optimal value of A is that which minimises the difference be- 
tween the functions E(N) for each sub-monolayer coverage 
in terms of least-squares. This method was not successful for 
all molecule-surface combinations, as there was not always 
a minimum in the least-squares calculation. In those cases 
where a minimum was calculated, the corresponding A value 
was found to be significantly lower than the expected value 
(typically around I0 6 s -1 , whereas the actual A expected 
was 10 12 ). This optimisation method for A assumes that 
there is no difference in the occupation of the high energy 
binding sites (i.e. the tail of the TPD spectra) with increas- 
ing coverage. Thus, all E(N) calculated for a given A value 
should coincide at low coverage; in our experiments, there 
were some minor differences in the tails of the TPD spectra, 
and consequently this method was rejected. To reproduce 
the TPD spectra at sub-monolayer coverages A was fixed at 
10 12 s- 1 . 

To successfully reproduce all the TPD data, from sub- 
monolayer to multilayer regimes, the results from the two 
coverage models are then combined, such that (a) when 9 ^ 
1, in Equation [5] E a d s is replaced by E(N) (derived from 
Equation [5}, n = 1, and N ^ N mono and (b) when 8 > 1, 
the TPD is the sum of the monolayer trace generated in (a), 
plus Equation^ where E a d s is taken from Table[5J A defined 
by Equation [3] n — and N is redefined as N — N mono . 

3.3 Discussion 

From Figure [3] it is clear that the combined model can very 
accurately reproduce the desorption profiles of all molecular 
species on all surfaces (O2, CO and CO2 on H2 0( np ), EbOjc) 
and SiO^). In particular, the model faithfully reproduces the 
sub-monolayer coverage TPD spectra. Although in a few 
cases the model data do not perfectly describe the empirical 
line profiles, the two clearly coincide within the 3 a errors on 
the model, as illustrated by the overplotted dot-dashed lines 
on the sub-monolayer coverage of CO2 on H20( c ). These 
errors are omitted from all the other plots for clarity. The 
errors result from the assumption that all the TPD spectra 
are concurrent in the highest energy binding sites, i.e. their 
tails overlap, which is a necessary assumption when applying 
this model. 

The multilayer TPD spectra, and their associated fits, 
are shown in the insets in Figure [3] Most importantly, as 



© 2002 RAS, MNRAS OOO.riTOI 



Thermal desorption characteristics of CO, 2 and C0 2 on H 2 0( np ), H 2 0^ and SiO x . 7 

Table 3. Calculated coefficients to describe the E(8), as in Equation [5] 



Molecule Surface Coefficients 









ai 


a 2 


a 3 


« 4 


c*5 


«(i 


a- 




«8 


a<) 


"10 






0.080 


0. 


0. 


0.030 


-4.124 


0.047 


-441.083 


-3.397 xlO 


-07 


8.449 


0.006 


-77.198 


o 2 


H 2 (c) 


0.099 


-0.032 


0.017 


0.043 


-306.959 


-0.098 


-13026.500 


-2.752 xlO 


-07 


8.212 


0.016 


-17.518 




SiO* 


0.028 


0. 


0. 


0.032 


-7.405 


0.027 


-118.135 


-3.323 xlO 


-09 


9.398 


0.067 


-0.248 




H20( np ) 


0.073 


0. 


0. 


0.054 


-428.915 


0.042 


-0.988 


-8.068 xlO 


-06 


7.447 


0.009 


-14.156 


CO 


H 2 (c) 


0.107 


-0.028 


0.011 


0.029 


-342.873 


0.010 


-26.879 


-4.150 xlO 


-07 


9.984 


0.010 


-0.585 




SiO* 


0.057 


0. 


0. 


0.060 


-1.176 


0.034 


-10.199 


-2.716 xlO 


-08 


9.571 


0. 


0. 




H20( np ) 


0.195 


-0.015 


0.007 


0.117 


-87.247 


-0.082 


-46.774 


5.035 xlO- 


07 


9.344 


0.037 


13.799 


co 2 


H 2 0( c ) 


0.201 


0. 


0. 


-32.300 


-6865.170 


0.026 


-73.854 


5.640 xlO" 


10 


15.280 


0.024 


-4.143 




SiOx 


0.199 


0. 


0. 


0.082 


-3.283 


-0.045 


-32.504 


-1.640 xlO 


-05 


5.937 


0. 


0. 



noted previously from the experimental data, the model data 
also show that multilayer desorption begins at the same tem- 
perature for each molecule-surface combination (the leading 
edges of the TPDs are identical) , although the onset of this 
desorption occurs at a higher temperatures on the R~2 0( c ) 
surface when compared with the amorphous H2 0( np ) and 
SiO x surfaces. From Table [2] it is clear that E a ds is consis- 
tently marginally greater on the crystalline surface, as com- 
pared to the amorphous surfaces, resulting in slightly lower 
desorption rates from crystalline surfaces, and the later on- 
set (in time or temperature) of the desorption process. This 
trend is even consistent with previously reported results (see 
Table [2]). Given that species such as CO, CO2 and O2 could 
be present on interstellar ices or grains at multilayer cover- 
ages of up to a few ML, this result indicates that a similar 
effect from underlying surface crystallinity could be relevant 
in interstellar regions. 

It should be noted that the desorption energies of 
species from H20( c ) are consistently higher than for the 
same species desorbing from either of the amorphous sur- 
faces, EbOfnp) and SiCv One possible explanation for this 
trend is that the ordered structure of the crystalline water 
surface induces larger sc ale interactions among the desorbed 
molecules. The study of |Dohnalek et al.ll 19991 ) revealed that 
amorphous water deposited on top of crystalline water dis- 
played a degree of crystalline nature in its first layer due to 
interaction with the ordered surface of the underlying crys- 
talline substrate. It is possible that O2, CO and CO2 experi- 
ence similar interactions with the crystalline water surface, 
although these bonds would likely be of lower energy than 
for a water- water system. Another possibility is that, rather 
than the crystalline surface causing unusual desorption en- 
ergies, it is in fact the values from amorphous surfaces that 
are altered due to the presence of defects in the amorphous 
structure which prevent normal monolayer and multilayer 
regimes to build up on these surfaces. 

Surprisingly, given its ubiquity in interstellar ices, the 
multilayer desorption energy of pure CO2 ice is not widely 
reported in the literature. One value has be en measured pre- 
vious ly from multilayers ofC02 onHOPG (jBurke fc Brownl 
2010), where analysis proved that the reaction order is zero, 
and estimated the magnitude of the pre-exponential factor 
at around 10 30 . However, consistent with the rather com- 
plex TPD spectra shown in grey in the insets on Figure [1] it 
was very difficult to undertake a detailed analysis of multi- 



layer CO2 desorption from th e HOPG surface. Conseq uently, 
the value of E a£ js quoted by iBurke fc Brown! (|2010T ) is not 
only somewhat higher than the value calculated here, but 
also has a quoted error greater than the value of E n d s itself 
(± 19484 K , which is omitted from Table [2] for clarity). 
Three other multilayer CO2 desorption measurements are 
claimed in the literature, as indicated in Table [2] unfortu- 
nately all three were determined from transmission infrared 
data, assuming the desorption kinetics were first and not 
zeroth order. Consequently, although included for complete- 
ness, as the values are often quoted in models, the numbers 
are not comparable to the results quoted here. 

For multilayer coverages of CO and O2 , the E a d s calcu- 
lated from this model can be compared directly to previous 
(multilayer) results in the literature (see Table [2}. Within 
the errors, our results for CO and O2 on H20( c ) - a crys- 
tal line surface - ar e ident ical to those obtained previously 
by lAcharvva et all (|2007m for 2 and CO on a poly crys- 
talline Au substrate, using the same fitting method. This 
last point is particularly r elevant when comparin g the re- 
sults, as both here and in lAcharvva et al.l l|2007h the pa- 
rameterisation of the fit is reduced to one by expressing the 
pre-exponential factor in terms of E a d a - Consequently, our 
pr e-exponential fac t or als o correlates well with that reported 
by lAcharvva et all |2007h previously. Even though the cal- 
culated E a ds for the amorphous surfaces are slightly lower 
than the result form the EbOfc) surface, they also are the 
consistent, with i n the calculated error, with the values from 
lAcharvva et al.l l|2007h . The other values of E a d s listed in Ta- 
ble [2] for CO and O2 on polycrystalline Au (ICollings et al. 
2003bl ; I6berg et all 120051 ; iFuchs et al.ll2006l ; iBisschop et~ 
2006T ). were calculated assuming a fixed pre-exponential fac- 
tor (A) of magnitude around 10 30 . As A and E a ds are not 
mutually exclusive, it is not possible to make a direct com- 
parison between these data and our results; however it is in- 
teresting to note that despite a difference of over four orders 
of magnitude in the value of the pre-exponential factor, the 
reported values of E a d s are all remarkably consistent. The 
question is only the extent to which these subtle differences 
in the value of TZ a ds affect the desorption of ice multilayers 
in interstellar space, and how to account for such subtleties 
in gas-grain models (see §2}. 

In modelling the sub-monolayer desorption of key inter- 
stellar molecules from different surfaces, the aim was to pro- 
vide the chemical modelling community with a useful equa- 
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Figure 3. Modelled temperature-programmed desorption spectra. Experimental data, as in Figure [TJ are plotted in black diamonds; 
coverages up to and including 1 ML are shown in the main window, with coverages of 1 ML and greater inset. Model data is overplotted 
as coloured/grey lines. Typical 3cr errors on the model are overplotted as dot-dashed lines on the plot of CO2 on H20( c ) only, for clarity. 
Some of the original TPD data are omitted for clarity (see Table [TJ. See text for full discussion. 



tion (Equation[5]and Table[3]) for mathematically describing 
the desorption characteristics of the surfaces, as well as an 
indication of the coverage at which a particular adsorbate's 
behaviour switches from the sub-monolayer to the multilayer 
regimes (a measure of N morlo ). The excellent agreement be- 
tween experimental data and the model is also validation of 
our method, described in § 13.21 The only limitation of this 
method is the certainty with which we have empirically se- 
lected the TPD curve used to model sub-monolayer desorp- 
tion characteristics. This curve represents the highest cover- 
age at which only monolayer desorption characteristics are 
seen, rather than, necessarily, a coverage of 10 15 molecules 
cm -2 (the generally accepted definition of one monolayer). 
However, it is possible that an exposure time of a few sec- 
onds more or less could have been a slightly better model 
system to represent (sub-)monolayer desorption. This is the 
limitation introduced by our chosen method of deposition 
calibration. 

Of all the data fitted here only the CO-ice systems 



show minute deviations at the highest sub-monolayer cover- 
ages, and such effects could be due to CO, as the smallest 
molecule studied, somehow being able to probe defects and 
cracks in the ice surface, prior to forming a m ultilayer, as de- 
scribe d previously in studies of CO on ASW jCollings et al.l 
l2003al lbfi, an effect omitted from this model. This may also 
explain why in these systems it is difficult to reproduce a 
perfect match in the transition regime between the mono- 
layer and multilayer coverages. Likewise, although the fits 
shown here to the C02-ice systems are excellent, a number 
of 'intermediate' TPD spectra, shown in grey in Figure [TJ 
with double and triple TPD peaks, had to be omitted from 
the modelling process. These spectra are indicative of 2D is- 
land growth on a surface, as seen in the formation of N2 ices, 
and indicates that the growth mechanis m of CO2 on any ice 
surface is different to that of CO or O2 (|6berg et al.ll2005l ; 
iFavolle et afll201ll) . Whilst each island is entirely indepen- 
dent of another on the surface and there is no diffusion be- 
tween them or bridges linking them, we can consider that we 
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Table 4. Calculated values of £(#) at a range of sub-monolayer 
surface coverage values, calculated using Equation [5] and the co- 
efficients from Table [3] 



Mol. 


Surface 


E(0) 1 K 














£(0.1) 


£(0.2) 


£(0.5) 


£(0.9) 


£(1.0) 






1161 


1082 


972 


928 


914 


o 2 


H 2 (c) 


1149 


1092 


1017 


975 


969 




SiO* 


1255 


1146 


1019 


945 


930 




H20( np ) 


1307 


1247 


1135 


956 


863 


CO 


H 2 (c) 


1330 


1288 


1199 


1086 


1009 




SiO* 


1418 


1257 


1045 


896 


867 




H20( np ) 


2346 


2258 


2197 


2197 


2236 


co 2 


H 2 (c) 


2514 


2451 


2364 


2341 


2361 




SiO* 


3008 


2798 


2487 


2317 


2271 



are in the sub-monolayer coverage regime. Once the islands 
have fully merged and overlayers of CO2 ice form, then we 
can be certain we are in a multilayer regime. This is what we 
assume here; again a more full description of the 'intermedi- 
ate' stages is beyond the scope of this paper. Nevertheless, 
such effects should not detract from the broad applicability 
of the sub-monolayer model. Furthermore, as a general rule 
the outcome of our monolayer TPD choice suggests that for 
CO and O2 it is reasonable to assume N mono = 10 15 , whereas 
for CO2 it is around N mono = 3 x 10 14 . This result is con- 
sistent with the exposures discussed previously in Figure [2] 
and reflects the CO2 island growth. 

As can be seen from the main panels in Figure [3] at 
very low coverages (< 0.5 ML) there is sometimes a small 
deviation in the leading edge of the modelled desorption 
compared to the experimental results, most evident on the 
SiOx surface. Although this lies within the error range of the 
fitting, there are two potential explanations for why such 
discrepancies could arise, both related to the assumption 
made in the model that all desorptions from a given surface 
share the same TPD trailing edge profile. As the sites in 
the tail of the TPD are the most energetically favourable, 
it follows that any adsorbate able to diffuse across the sur- 
face will bind most readily to these sites, occupying them 
first on adsorption and last on desorption. Such an effect 
has clear ly been seen previo usly with H2 and D2 on ice 
surfaces (|Amiaud et al J 12006). Even under ultra-high vac- 
uum conditions, molecules such as H2, CO and H2O are still 
present in the experimental chamber, all be it at exception- 
ally low concentrations (partial pressures around 10~ 13-14 ). 
Consequently, over the very long timescales required to com- 
plete the experiments described here (14-16 hours per single 
molecule-surface combination), it is potentially possible that 
some of the surface sites could be dynamically occupied for 
a short period by a 'pollutant' adsorbate, thus blocking the 
adsorbate molecule of interest from occupying a highest en- 
ergy binding site. Such an effect would then lead to higher 
occupation of the lower energy binding sites and a miss 
match between the experimental and model data, the lat- 
ter appearing at slightly higher temperatures. Such an effect 
would not have been evident in the work of lAmiaud et al.l 
(2006) because they were only cyclically heating their ice 
surfaces between 10 and 30 K, and consequently could com- 
plete the experiment on a much shorter timescale than was 



possible here. As H2 and CO are reasonably volatile pollu- 
tants, the build up of H2O would be most likely to have an 
effect on the nature of the surface, and such an effect would 
be most evident on the non ice surface, i.e. SiO^. However 
as the experimental sequence was always run from low to 
high coverage, if this were the key reason for any discrep- 
ancy, it should be more pronounced as we get towards the 
monolayer coverage data. Evidently this is not the case. 

It is also important to consider mass effe cts; in compar- 
ison to the H2-ice system studied previously (|Amiaud et al.l 
2006), CO, O2 and CO2 are more massive molecules, and 
therefore likely to diffuse more slowly on the surfaces. There- 
fore, unlike the H2, when deposited on the surface at rela- 
tively low concentrations, CO, O2 and CO2 may not sam- 
ple all the binding sites before adsorbing, just ballistically 
depositing in a 'stick and stop' process. Subsequently, dur- 
ing the TPD warming, all the molecules will diffuse, but 
again the heavier molecules may not fully sample the surface 
before desorbing. Then, at very low coverages, not all the 
highest energy binding sites will be occupied, and the TPD 
tails may not be quite coincident. As a result, the model 
would slightly underestimate the leading edge of such TPD 
curves. From this explanation it is also possible to rationalise 
why the discrepancy between the model and empirical data 
would be greatest on the SiOa; surface: the breadth of the 
sub-monolayer TPD peaks from the molecule-silicate sys- 
tems in comparison to the molecule-ice systems clearly indi- 
cates a broader range of binding sites on the SiO^ surface. By 
contrast, the water surfaces have a narrow, and an energeti- 
cally similar range of binding sites, which may also intimate 
that the ice surfaces do not ha ve many dangling bonds or 
much proton disorder in them |Fraser et alj 12004 ). This is 
contrary to what we might expect, as it is often assumed that 
amorphous water ice potentially has a very broad range of 
binding sites at its surface. However without spectroscopic 
data we cannot comment further on this here. Consequently, 
the probability of sampling the highest energy binding sites 
is even lower at low coverages on SiOa, than the ice surfaces, 
provided the diffusion rates on all three surfaces are approxi- 
mately equal. Nevertheless, these tiny differences should not 
detract from the excellence of the fits, nor that the model 
can be executed whilst ignoring the pumping speed in the 
chamber; this is because in FORMOLISM the pumping is so 
effective that as soon as a molecule desorbs from the surface 
it is effectively removed from the chamber - in experiments 
without such efficient pumping the experiments described 
here would simply be impossible to undertake, as the TPD 
tails may never overlap. 

To compare the E{ff) calculated here with previous ad- 
sorption energy values of CO, O2 and CO2 as reported in the 
literature, we have produced Tabled This shows the specific 
values of E{6) for each molecule on each surface at coverages 
ranging from 9 — 0.1-1.0. The errors on each energy value 
are determined by the error on the fit to the experimental 
data and are around ± 25 K in all cases. The comparison 
is rather complicated by the fact that no-one else has ever 
considered fitting coverage dependent adsorption energies 
for these particular molecules at sub-monolayer coverages 
on surfaces of astrophysical relevance. In general the effec- 
tive adsorption energy is higher at lower coverage, reflecting 
that the highest energy binding sites are occupied first, and 
tends towards the calculated multilayer value for the same 
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molecule-surface system as the coverage rises towards that 
of a monolayer. Although, by the monolayer coverages, the 
trend again emerges that the binding energy on the crys- 
talline surface is consistently higher than on the amorphous 
surfaces; the same is not true at lower coverages. 

Previous studies using different methods have re- 
ported t he monolayer energy of CO as: 1564 ± 120 K on 
HOPG (luibricht et al.ll2006T ). 1624 ± 360 K on meteorite 
(iMautner et al.ll2006h. and 1179 ± 24 K on a highly porous 
ASW (|Collings et al.ll2003b| y Certainly our values for CO on 
the ice surface s cover a range of values including the result of 
ICollings et all j2003rJ ). and at the lowest coverages our value 
on the Si02 surface is closer to that reported for HOPG. it 
is pertinent to compare the lowest coverage value since the 
HOPG data were measured in a molecular beam scatter- 
ing experiment, indicating that only a low concentration, 
dynamic coverage was resident on the HOPG surface for a 
relatively short time, and therefore must have been occupy- 
ing the highest energy binding sites. The meteorite sample 
would no-doubt encompass a much wider variety of binding 
sites than event the amorphous silicate used here, so again 
it is not surprising that the desorption energy is even higher 
on this surface. This excellent agreement between our and 
previous results is further justification that our model is vi- 
able and effectively reproduces the effects of sub-monolayer 
desorption from a variety of interstellar-relevant surfaces. 

For O2, previous calcu lated sub-monolayer energies are: 
1082 ± 1 20 K on HOPG (lUlbricht et al.ll2006T ). 1203 K on 
graphite (|Boian fc Steele! I1987T ). Again our values at the 
lowest coverages on Si02 compare favourably with previ- 
ously reported values from carbon-based surfaces. The des- 
orption energy of the CO 2 monolayer has been calculated as: 
2766 ± 241 K on HOPG JUlbricht et al.ll200fj) . 3079 ± 20 K 
on grap hite (iTerlain fc Larherl Il983l ') 2553 ± 232 K on 
H 2 (c) (|Andersson et al.ll2004l), 2490 ± 240 K on H 2 (np) 
and 2393 ± 240 K on H^Or.) (|Galvez et al.ll2007h . Certainly 
the range of our data compare favourably with these values; 
in particular at the lowest coverages, the value of E{9) for 
CO2 on Si02 is consistent with the molecular beam scatter- 
ing result from HOPG, and the value of E{6) for CO2 on 
H2 0f e i is identical within t he quoted errors to that measured 
bv lAndersson et al.l (|2004r ). who also used a molecular beam 
scattering experiment and therefore probed the highest en- 
ergy binding sites for CO2 on the cr ystalline ice surface . In- 
terestingly, the value calculated by (jGalvez et al. I [20071 ) for 
desorption from H20( np ) is higher than that from H20( c ), 
a result in contrast to the data presented here. However, 
within the stated experimental errors of the study, these 
values are consistent with our results. Such excellent over- 
lap between the results also verifies our earlier assumptions 
in terms of choosing the appropriate TPD spectrum to de- 
fine as the monolayer coverage for CO2 on the ice surfaces, 
and indeed our calculation of N morlo . 

The adsorption behaviour of the molecule-surface com- 
binations can be divided into roughly three categories re- 
lated to the difference between the multilayer energy (or 
adsorbate-adsorbate interactions) and the energy distribu- 
tion of the surface adsorption sites. These behaviours can be 
explained with reference to Figure [4] a plot of the modelled 
energy as a function of coverage, E(8). Where the energy of 
the multilayer is significantly higher than the energy distri- 
bution at coverages below 1 ML, as is the case for CO on 
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Figure 4. Total modelled energy as a function of 
exposurc/N mo ,j for all molecule-surface combinations. Plots 
are as follows: O2 dot-dashed, CO solid, CO2 dashed; H20( np j 
black, H20( c j dark grey, SiOx light grey. 



H20( c ), molecules will preferentially adsorb onto the sub- 
strate at low coverages, rather than starting to form a mul- 
tilayer or island. In this 'full wetting' regime, the multilayer 
appears only after all of the sites on the surface are filled. 
All the CO lines (solid) in Figure [4] follow this regime. 

Where the multilayer energy is lower than the mono- 
layer energy, such as for CO2 on H20( np ), molecules prefer- 
entially form clusters (or islands) on the surface at very low 
coverages, well before 1 ML. In this 'non-wetting' regime, 
these islands, which are multilayers thick, do not cover the 
whole substrate but eventually build long towers which may 
join and overlap, until eventually multilayer zeroth order 
desorption is observed. All the CO2 sub-monolayer cover- 
ages on any H2O ice surface follow this regime, see Figure [4] 

These two 'wetting' and 'non-wetting' regimes can be 
further contrasted by considering the TPD spectra of CO 
on H20( c ) versus those of CO2 on H20( np ) in Figure[T] CO2 
rapidly forms a double-peaked multilayer desorption profile 
while CO shows a gradual filling of the surface sites before 
a multilayer feature forms. 

Intermediate between these two regimes, molecules first 
fill the lowest energy surface sites, before starting to form 
multilayer islands on top of the existing monolayer, before 
the monolayer coverage is fully complete. This occurs where 
the energy of the final binding sites in the monolayer is equal 
to or greater than that of the multilayer binding sites, and 
usually is only seen at exposures close to N mono . These three 
states are relatively easy to observe in Figure[4] although the 
exact boundaries between the regimes are difficult to define, 
as island growth is difficult to reproduce and analyse ex- 
perimentally. Nevertheless, the regimes associated with each 
molecule-surface combination are summarised in Table [5] 
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Table 5. Adsorption behaviour of each molecule-surface combi- 
nation, estimated from the modelled energies in Tables [2] &|3] and 
the desorption characteristics in Figure [l] a 



Molecule 


Surface 


Behaviour 




H20( np ) 


Full wetting / Intermediate 


2 


H 2 (c) 


Full wetting / Intermediate 




SiO* 


Full wetting 




H20( np ) 


Full wetting 


CO 


H 2 {c) 


Full wetting 




SiO x 


Full wetting 




H20( np ) 


Non-wetting 


CO2 


H 2 (c) 


Non-wetting 




SiOx 


Non-wetting / Intermediate 



a 'Full wetting' behaviour is seen for molecules which fill all the 
lowest energy sites on the surface before forming a multilayer. 
'Intermediate' behaviour is when the molecules start to fill the 
surface, but form islands before it is full. 'Non-wetting' behaviour 
is when islands are almost immediately formed by molecules on 
the surface. 



Temperature / K 

30 40 
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Figure 5. The simulated desorption of 2 ML of CO (solid), 
O2 (dot-dashed) and CO2 (dashed) from H 2 0( np ) (black), 
H20( c ) (dark grey) and SiO^ (light grey) at a heating rate of 
1 Kcentury -1 . For comparison, the black dotted line shows the 
modelled desorption characteristics of CO on H20( np ) if only 
multilayer desorption is included in the model. 



4 ASTROPHYSICAL IMPLICATIONS 

Under astro-physical conditions, the heating rate is likely 
to be on the order of 1 Kcentury -1 (hot core heating rate 
jViti fe Williams! 1 19991 )1. many orders of magnitude slower 
than the laboratory value of 10 Kmin - employed here. By 
simulating the desorption profiles of the molecule-surface 
combinations here, it is possible to investigate whether sur- 
face type and coverage is likely to affect desorption charac- 
teristics on astrophysically relevant timescales. 

Figure [5] presents the results of such a simulation for 
all molecule-surface combinations, assuming a 2 ML cover- 
age of each species on each grain (i.e. one multilayer and 
one monolayer). This seems reasonable as none of these 
species (except possibly CO) form large pure ice overlay- 
ers in the ISM. In Figure [SJ the dot-dashed line shows the 
modelled desorption characteristics of 2 ML of CO from 
H2 0( n p) if only multilayer desorption is considered. It is ev- 



ident that without including the sub-monolayer component 
of the model, the desorption time is vastly underestimated 
(by over a few x 10 4 years or 10 - 15 K), particularly for 
CO and CO2 . This is a comparable error to that introduced 
when gas-grain models treat m ultilayer desorption as first 
order, as suggested previously (|Fraser et al.|[200ll ). and il- 
lustrates that sub-monolayer desorption characteristics are 
fundamental to describing the overall desorption profile of 
species on astrophysical timescales. 

Figure shows that for O2 there is the least differ- 
ence between the time (or temperature) at which desorption 
would be completed in the multilayer only model versus the 
combined multilayer-monolayer model. However, monolay- 
ers of both CO and CO2 are able to reside directly at the 
surface for thousands of years, or at equivalent temperatures 
10 - 15 K higher, than those predicted by multilayer desorp- 
tion alone, particularly on the amorphous silicate surface. In 
quiescent regions of the interstellar medium or at the edges 
of molecular clouds, where grains are warmer, this would 
suggest that small molecules, like CO and CO2, could eas- 
ily be present on the grain surfaces, potentially undergoing 
chemistry to generate larger chemical species, long before 
detectable ice layers r eally form. Such p rocesses have been 
postulated previously (|Fraser et al.1120051 ). and it remains an 
interesting follow-on from this work to see which chemical 
species will form in astrochemical models when these sub- 
monolayer desorption data are considered. 

What is also clear from Figure [S] is that, while the 
molecules desorb at different times (and temperatures), 
there is a clear pattern to the order of desorption from 
the various substrates. The onset of desorption from the 
amorphous surfaces, H2 0( np ) and SiOa;, occurs almost con- 
currently for all species, while desorption from HbOfc) oc- 
curs later, at a higher temperature. This delay in des- 
orption from crystalline surfaces indicates that in regions 
where crystalline water is predicted to domina te, for exam- 
ple o n processed ices inside disks aro und YSOs dVisser et all 
2009), or even cometary surfaces iRoushl [iooii r volatile 
species will not start desorbing from grain surfaces for 
a few thousand years longer (or over temperature ranges 
of a few kelvin higher) than anticipated from zeroth or- 
der desorption kinetics models. Such differences are sub- 
tle, but when one considers that these molecules will be 
highly mobile at such temperatures, the potential for com- 
plex chemistry to occur on the grain surfaces is increased. 
Furthermore, it is exactly these inner disk regions, where 
crystalline water ice is predicted to be found, from which 
hot core and hot corino gases are populated by desorb- 
ing molecules(e.g. IViti fc Williamslll999l; [van der Tak et al 



pOtt ; ICazaux et al.ll2003l ; ISemenov et al.l l2010l ; iFavre et al 
hoill ). 



Similarly, the simulation results in Figure [S] show the 
residence time of molecules on a H20( c ) surface is signifi- 
cantly higher than on the amorphous surfaces. Starting from 
a 2 ML coverage, the time taken for 50 % of the O2 on SiO^ 
or H20( np ) to desorb is 8000 years, while on H20( c ) this 
value is 8700 years, an increase of 8.75 %. Similarly, the dif- 
ference between desorption of 50 % of the molecules from 
H 2 {np) and from H 2 (c) is 8.82 % for CO and 6.97 % for 
C0 2 . 
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5 CONCLUSIONS 

A key aim of this work was to address what effect the un- 
derlying dust grain surface could have on the desorption 
characteristics of molecules adsorbed there, and how this 
desorption is modified in the sub-monolayer coverage regime 
compared to the multilayer regime. 

This work used temperature programmed desorption 
to measure the desorption characteristics of O2, CO and 
CO2 from H20(, lp ), EbOjc) and SiOx over sub-monolayer 
to multilayer coverages. The experimental data were mod- 
elled using the Polanyi-Wigner equation, combining differ- 
ent approaches to reproduce the sub-monolayer and multi- 
layer TPD traces, in particular by varying F, a ds as a func- 
tion of coverage in the sub-monolayer regime. The desorp- 
tion can be categorised as 'full wetting', 'intermediate' or 
'non-wetting' behaviour, and the switching point is well de- 
scribed by N mono , an empirical measure of the coverage at 
which pure monolayer desorption stops. On both a labora- 
tory and an astrophysically relevant timescale, the desorp- 
tion characteristics of molecules from the amorphous sub- 
strates H2 0( np ) and SiO x were found to be very similar, 
while on the crystalline surface H2 0( c ), molecules desorbed 
at higher temperatures and on a longer timescale. 

It would be relatively trivial to implement the data pre- 
sented here into astrochemical models; it would require only 
a subroutine to describe the sub-monolayer, coverage depen- 
dent, E ac j s , which can be calculated using the parameters 
presented in this paper, and the coverage N mono , at which 
the true monolayer desorption ends and beyond which a 
fixed E at j s , and zeroth order desorpt ion can be implemented. 

Previously. iGreen et a l. (2009) concluded that surface 
type is relatively unimportant, and rather it is the heating 
rate and grain size that dominate desorption kinetics; we do 
not disagree that both these factors are important, but by 
not considering the sub-monolayer case, this previous work 
overlooks certain intricacies. Firstly, the results presented 
here suggest that the surface type is the dominant factor 
controlling desorption under sub-monolayer coverage condi- 
tions. This definition of surface type must include both the 
surface material and the degree of crystallinity, and not sim- 
ply the material alone. From the data presented here, it is 
clear that the desorption characteristics of molecular species 
from H20( np ) and SiOx, two amorphous surfaces, are much 
more similar than those from H20( np ) and H20( c ), which are 
composed of the same underlying material. Secondly, while 
the size of the grain is critical to desorption characteristics, 
it is both the surface area and the surface-adsorbate interac- 
tion which must be considered. Without including informa- 
tion on whether the surface-adsorbate interaction is wetting 
or non-wetting, it is not possible to know the coverage at 
which behaviour switches from sub-monolayer (to interme- 
diate) to the multilayer regime, i.e. the value of N morlo . As 
it is notoriously difficult to calculate physisorption surface- 
adsorbate energies theoretically, this is necessarily an em- 
pirically determined value, highlighting the importance of 
experiments such as those presented here. 
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